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Abstract Cosmic ray acceleration in SNRs in the presence of the Alfvenic drift is considered. It is 
shown that spectra of accelerated particles may be considerably softer in the presence of amplified 

magnetic fields. 

Introduction 

It is almost doubtless now that supernova remnants (SNRs) are the main source of galactic cosmic 
rays (CR). The outer shell of the exploding star moves with a supersonic velocity and produces a strong 
shock wave in the circumstellar medium. Diffusive shock acceleration [TJ [5] results in the energy gain of 
i^ ' energetic particles. The observation of high-energy TeV gamma-rays from several SNRs is the evidence 

of effective acceleration of cosmic ray particles up to the energy about 100 TeV [3] in these objects. This 
energy is only a factor of 30 lower than the "knee" energy Eknee ~3 PeV in the CR spectrum. 

At present there exist two numerical models of nonlinear diffusive shock acceleration at the moving 
spherical shock of Berezhko and co-authors [1], and Jones and co-authors [S]. A so-called shock mod- 
ification by the CR pressure is taken into account in these models. This is important since about 10 
percent of supernova energy is transferred to accelerated CRs, if SNRs are the main source of CRs in the 
Galaxy. In addition it seems that CR particles are accelerated only at some part of the SNR shock, as it 
1-^ ' is observed in several SNRs. 

'j^l These models may be used for modeling of CR acceleration in particular SNRs and for the calculation 

Q , of an overall CR spectrum that is produced during lifetime of a SNR. It is expected that this spectrum 

$-H ■ is close to _E^^ (see j^). Slightly harder spectrum E~^-^ was numerically obtained recently [7j. 

lyj I In this short report we present results of the new numerical model of CR acceleration in SNRs. It 

C^ . is shown, that CR spectra may be significantly softer if the advective velocity of CRs downstream of 

the shock is essentially different from the gas velocity. This situation is probable if the magnetic field is 
amplified in SNRs. 

'.^ ' Model of nonlinear diffusive shock acceleration in SNRs. 

l/~j I Hydrodynamical equations for the gas density p(r, t), gas velocity u(r, t), gas pressure Pg{r,t), and 

[~^ . the equation for isotropic part of the CR momentum distribution N{r, t,p) in the spherically symmetric 

Ps| ' case are given by 
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Here Pc = An J p'^dpvpN/3 is the CR pressure, w{r, t) is the advective velocity of CRs, ^g is the adiabatic 
index of the gas, and D(r,t,p) is the CR diffusion coefficient. It was assumed that diffusive streaming 
of CRs results in the generation of magnetohydrodynamic (MHD) waves. CR particles are scattered by 
these waves. That is why the CR advective velocity w may differ from the gas velocity u. Damping of 
these waves results in additional gas heating. It is described by the last term in Eq. (3). The last term 
in Eq. (4) corresponds to the injection of thermal protons with momentum p = pinj and mass m at the 
shock front at r = R(t). The dimensionless parameter rj determines the injection efficiency. 

CR diffusion is determined by magnetic inhomogeneities. Strong streaming of accelerated particles 
changes medium properties in the shock vicinity. CR streaming instability results in the high level of 
MHD turbulence [2] and even in the amplification of magnetic field in young SNRs [8 . Due to this effect 
the maximum energy of accelerated particles may be higher in comparison with previous estimates [9] . 

According to the recent numerical modeling of this instability, magnetic field is amplified by the flux 
of run-away highest energy particles in the relatively broad region upstream of the shock \Wl ■ Magnetic 
energy density is a small fraction {^ 10"'^) of the energy density of accelerated particles. This amplified 
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Fig.l. Dependencies on time of the shock 
radius R (thick sohd hne), shock velocity R 
(thick dashed hne), the total compression ra- 
tio of the shock a (dotted line). Dependencies 
of the magnetic field strength downstream of 
the shock (dashed line) and ratio of CR energy 
and energy of supernova explosion Ecr/EsN 
100 (solid line) are also shown. 



almost isotropic magnetic field can be considered as a large-scale magnetic field for lower energy parti- 
cles which are concentrated in the narrow region upstream of the shock. Streaming instability of these 
particles produces MHD waves propagating in the direction opposite to the CR gradient. This gradient 
is negative upstream of the shock and MHD waves propagate in the positive direction. The situation 
changes downstream of the shock where CR gradient is as a rule positive and MHD waves propagate in 
the negative direction. This effect is mostly pronounced downstream of the shock because the magnetic 
field is additionally amplified by the shock compression and the Alfven velocity Va = B/y/Airp may be 
comparable with the gas velocity in the shock frame u' = R— u{R— 0, t). As for CR diffusion coefficient, 
it is probably close the Bohm value Db = pvc/3qB, where q is the electric charge of particles. 

Numerical modeling of CR acceleration in SNRs. We apply finite-difference method to solve 
Eqs (1-4) numerically upstream and downstream of the shock. The auto-model variable i^ = r/R{t) is 
used instead of radius r. The non-uniform numerical grid upstream of the shock at r > R allows to 
resolve small scales of hydrodynamical quantities appearing due to the pressure gradient of low-energy 
CRs. Eq. (4) for CRs was solved using an implicit finite-difference scheme. The explicit conservative 
TVD scheme |11) for hydrodynamical equations (1-3) was used. These equations are solved upstream the 
shock using the explicit finite-difference scheme. 

We shall assume that the coordinate dependencies of the magnetic field and the gas density coincide: 



B = v/47rpo 
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(5) 



Here po is the gas density of the circumstellar medium. The parameter Ma determines the value of the 
amplified magnetic field strength. The magnetic energy is about 3.5 percent of the dynamical pressure 
PqR^, according to estimates from the width of X-ray filaments in young SNRs [T^]. This number and 
characteristic compression ratio of a modified SNR shock cr — 6 correspond to Ma ~ 23. 

CR advective velocity differs from the gas velocity on the value of the radial component of the Alfven 
velocity Vat calculated in the isotropic random magnetic field: w = uzL Vat- Here signs ± correspond to 
regions upstream and downstream of the shock respectively. Using Eq. (5) we obtain 
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We shall use CR diffusion coefficient D — Db calculated with the magnetic field strength (5). Though 
in the real situation the level of MHD turbulence may drop with distance upstream of the shock and 
diffusion may be faster than the Bohm one (see [ID]), we shall use these assumptions here (see also [7]). 

The numerical results are obtained for the SNR shock propagating in the medium with the number 
density uq — 0.1 cm~^ and the temperature T = 10* K. We use the ejecta mass M^j — 1.4M0, the energy 
of explosion Esn — 10^^ erg and the parameter of ejecta velocity distribution fc = 7, corresponding to 
la supernovae. The initial shock velocity is Vq =31000 km/s. The injection efficiency is taken in the 
form r] — 0.01i?/Vo, and the injection momentum is pinj — 2m{R — u{R 4- 0, i)). This dependence of 
the injection efficiency on the shock velocity results in the significant shock modification already at early 
stages of SNR expansion. This is in agreement with the observations of young extragalactic SNRs [13] 
and with the modeling of coUisionless shocks [Tl] . 

The numerical results are obtained using the uniform grid with 800 cells downstream of the shock for 
the variable ^, uniform grid with 800 cells upstream of the shock for the variable ln(<^ — 1 + 10^^^). The 
uniform grid with 200 cells for the variable \np/mc is used. 




Fig. 2. Radial dependencies of the gas den- 
sity (thick sohd hne), the gas velocity (dot- 
ted hne), CR pressure (thick dashed hne), the 



gas pressure (dashed line) at t 



10^ 



yr- 



At 



this moment of time the shock velocity is 3890 
km/s, its radius is 7.6 pc, the magnetic field 
strength downstream of the shock is 136 /iG. 
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Fig. 3. Spectrum of CR protons at the shock 
front at 10"^ yr after explosion. The results 
corresponding to the Eq. (6) (solid line) and 
to CR advection velocity coinciding with the 
gas velocity downstream of the shock (dashed 
line) are shown. 



The dependencies on time of the shock radius R, the shock velocity R, the total compression ratio of 
the shock a, the magnetic field strength downstream of the shock and CR energy Ecr/EsN are shown 
in Fig.l. The calculations are performed until the beginning of the radiative phase of SNR expansion at 
t — 10^ yr, when the value of the shock velocity drops down to i? = 206 km/s. At this moment of time 
the maximum energy of particles accelerated in SNR is about 10 TeV, while higher energy particles have 
already leaved the remnant. This maximum energy may be significantly lower if one takes into account 
the wave damping on neutrals or nonlinear damping |15j . 

Radial dependencies of physical quantities at t — 10'^ yr are shown in Fig. 2. The contact discontinuity 
between the ejecta and the interstellar gas is at r = 6.2 pc. The reverse shock in the ejecta is situated 
at r = 5.2 pc. We neglect the injection of thermal ions into diffusive shock acceleration at the reverse 
shock. At the Sedov stage the reverse shock moves in the negative direction and reach the center. The 
appearing reflected shock wave moves in the positive direction and overtakes the main shock wave after 
40 thousand years after explosion. This results in the non-monotonic behavior of the shock velocity at 
this moment of time in Fig.l. 

The spectrum of CR protons at the shock at t = 10'^ yr is shown in Fig. 3 (solid line). It is compared 
with the proton spectrum obtained when the CR advection velocity coincides with the gas velocity 
downstream of the shock (dashed line). CR pressure is Pc = 0.25pqR^ in the first case and a factor of 
two higher in the second one. The CR spectrum is significantly softer in the first case. 

The spectrum of CR protons produced during the lifetime of the SNR F{p) is shown in Fig. (4) (solid 
line). It is determined by the integration of CR momentum distribution N{p) on the simulation volume 
at the moment t = 10^ yr. It is compared with the CR proton spectrum obtained when the CR advection 
velocity coincides with the gas velocity downstream of the shock (dashed line) . The last curve is similar to 
results of Berezhko and Volk [T| . The integrated CR spectrum is again significantly softer in the first case. 

Conclusion. Magnetic field amplification in SNRs may result in the significant difference of the 
CR advective velocity and the gas velocity downstream of the SNR shock. CR spectra in young SNRs 
and CR spectra produced by SNR during its lifetime may be significantly softer due to this effect. The 
fraction of supernova energy transferred to CRs is also reduced from 75 percent down to 35 percent of 
the mechanical energy of explosion. SNR luminosity in the TeV gamma-rays is reduced also. 
This work was supported by RFBR 07-02-00028 grant. 
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Fig. 4. Spectrum of CR protons produced dur- 
ing the lifetime of SNR. The results corre- 
sponding to the Eq. (6) (solid Hue) and to CR 
advection velocity coinciding with the gas ve- 
locity downstream of the shock (dashed line) 
are shown. 
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